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ABSTRACT 



The low degree splittings obtained from one year 
of GOLF data analysis are combined with the MDI 
medium-1 144-day splittings in order to infer the so- 
lar internal rotation as a function of the radius down 
to 0.2i?Q. Several inverse methods are applied to the 
same data and the uncertainties on the solution as 
well as the resolution reachable are discussed. The 
results are compared with the one obtained from the 
low degree splittings estimated from GONG network. 
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1. INTRODUCTION 



Figure 1. Sideral sectoral splittings (GOLF) and ai coef- 
ficients (MDI and GONG) with their formal errors as a 
function of the turning points of the modes. 



The rotation of the solar core is an important open 
question that can be addressed by using spatial data 
from SOHO. In particular, GOLF experiment is ded- 
icated to the observation of low-degree oscillations 

which sound the core^ Here we use the GOLF 

frequ ency splittings ( poster 1.39, these proceedings 



19981) obtained from one year of observation begin- 



nmg on April 11th 1996 together with t he MDI 144- 
day s plittings for degree up to ^ = 250 (fechou ct al 
1998 ) . For comparison we have also used th e (t()N(i 
splittings of low-degree mo des obtained ( RabcUo- 



Soares fc Appourchaux 1998|) from 1 year of ground- 
based observations (August 1995- August 1996). 

Figure ^ shows the GOLF sectoral splittings with 
their formal errors whereas for MDI and GONG data 
the figure shows the n.i coefficients of the expansion of 



the splittings on [RitzwoUer fc Lavely (1991)| polyno- 



mials. These two quantities (sectoral splittings and 
Oi coefficients) may differ slightly in theory for I > 1 
because of the latitudinal dependence of the solar ro- 
tation. In Section |] we briefly recall how the 2D in- 
verse problem related to the internal rotation can be 
reduced to a ID problem for either the ai coefficients 
or the sectoral splittings and how in both cases the 
use of sectoral splittings may need some ajiriori as- 
sumptions on the rotation. Then, Section pj presents 
the different inverse methods used a nd, nnally, we 
discuss the results obtained in Section pi 



FROM THE 2D TO THE ID INVERSE 
PROBLEM 



For low rotation the frequency of a mode of radial 
order n and degree I is splitted in 2Z -I- 1 components 
of azimuthal order m = —I, .. -\- 1 and the splittings 



Im 



are given by: 



Knimir, fi)n{r, fi)drdn, (1) 



where n{r, /j) is the unknown rotation rate versus 
depth r and colatitude 9 {fi = cos 9) and Knim{r,^) 
the so-called rotational kernels calculated for each 
mode from oscillation eigenfunctions of an equilib- 
rium solar model. 

In order to investigate the rotation below 0.4i?Q 
where the latitudinal dependence is particularly not 
well constrained because of the few azimuthal orders 
provided by the low I degree modes, one may want to 
simplify the problem and reduce it to ID problem in 
radius. As a matter of fact 2D inversions are usually 
not able to peak averaging kernels in both radial and 
latitudinal direction below OARr^ and two different 
ID approximations of Equation ^ are used instead. 
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2.1. ID relation for sectoral splittings 

One possibility in order to obtain a ID integral re- 
lation is to search for the rotation rate at a given 
latitude only. Since the most constrained zone is the 
equatorial one, we try to investigate only the equa- 
torial rotation rate. For this purpose we use the ap- 
proximation 



(2) 



Cuypers 1980| ) 
Drmalized sucn 



where Kni{r) are radial kernels (e.g. 
and P["{fi) are Legendre functions normalized 

that J^-^ P[^{iiydii = 1, and which satisfy the fol- 
lowing property for m = I: 



{2l + l)\ 

22/+l(^!)2 



(3) 



This shows that, for high-degrees ^, the major con- 
tribution to sectoral splittings Ai/„n comes from the 

equatorial rotation rate f2ei3('') = M = !)• This 
leads to the ID integral approximation: 



The radial kernel Kni (r) is the same as in Equation || 
but, from Equation]^, the function f2i(r) obtained by 
inverting ai coefficients corresponds to the searched 
rotation rate only where the rotation do not depend 
on the latitude. Otherwise it corresponds to some av- 
erage over latitudes that can be estimated by looking 
at the corresponding 2D averaging kernel (cf. Sec- 
tion B] and Figure 0). 



3. INVERSE METHODS 



We have used two kinds of inverse methods for solving 
the ID integral equations. 



1. A Regularized Least-Squares (R.L S) method wit h 
Tikhonov regularization (see e.g. [Corbard 1997|) . 
This is a global method which gives a solution at 
all depths which fits the data at the best in the 
least square sense. This is a linear method and 
then the value of the rotation obtained at any 
radius ro is a linear combination of the data: 



Kni[r)Q.eq{r)dr 



(4) 



Nevertheless this approximation is valid only for high 
degrees I. For lower degrees the sectoral splittings 
are sensitive not only to the equatorial rotation but 
also to the rotation rate in a large angular domain 
around the equator. The extent of this domain and 
the influence of this approximati on on the estimation 
are discussed in Corbard (1997)1. 



2.2. ID relation for ai coefficients 

Some experiments like GONG and MDI produce a 
small number of the so-called a-coefficicnts of split- 
tings expansions on a set o f orthogonal polynomials 
(Ritzwoller & Lavely 1991). Assuming that the re- 
lation between individual splittings and these coeffi- 
cients is linear, an equation similar to Equation |l| can 
be established by computing the appropriate kernels 
K^i- (r, y) related to each a"'-c oefficients for odd in- 
dices j (see e.g [Piipers 1997). Furthermo re it has 



been shown by Ritzwoller & Lavely (1991) that the 



expansion of the splittings in orthogonal polynomials 
corresponds to an expansion of r2(r, /i) such that: 



j=3,5.. 



(5) 



y(^o) = ^Cni{ro) 

nl 



(7) 




By replacing in Equation ^ or ^ we obtain: 



^1 - ('^o) = / {}] Cni{ro)K^i{r) ] 17, {r)dr 

(8) 

The function in parenthesis is called ID averag- 
ing kernel at tq, ^(r, tq). The result obtained at 
ro will be easier to interpret as a local average of 
the rotation when this kernel is well-peaked and 
without strong oscillatory behaviour. 

2. Two 'local' methods which search directly the 
coefficients C„i(ro) which are able to peak the 
averaging kernel near tq. The two methods dif- 
fer essentially in the way to localize the averag- 
ing kernel. The SOLA (Sub tractive Optimally 



Localized Average) method ( |Pijpers fc Thomp 



son 1992 ) fits the averaging kernel to a Gaus 
sian iunction of given width w hereas the MOLA 



meth od (Multiplicative OLA) ( Backus fc Gilbert 
19701 ) simply gives high weights m the mmimiza- 



tion process to the part of the averaging kernel 
which are far from the target radius. In both 
case we use a regularizing parameter in order to 
establish a balance between the resolution and 
the error magnification reached at the target ro- 



where Pj are the Legendre polynomials. This forms 
the so called 1.5D problem where each Oj-coefficient 
is related to the expansion function of the same index 
through a ID integral. Therefore the first term of the 
expansion Equation H which do not depend on the 
latitude can be related to the ai-coefficients through: 



Kr,i(r)hi{r)dr 



(6) 



HOW TO USE / = 2, 3 SECTORAL 
SPLITTINGS? 



As already quoted, GOLF data for ^ = 2, 3 are not 
oi coefficients but sectoral splittings. Therefore one 
may want to use the Equation ^ in order to infer the 
equatorial rotation rate. There are two difficulties 
with this approach: 
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Ai.„„[n(r,e)=D,^(r)]-4F„„[Q(r,0)] 




^' """"" iiiiiiiiiiiiiiiiiiiiilliiiiiiiiiiiiiiiiiiiillii- 
200 250 



Figure 2. Difference tn sectoral splittings computed from 
Equation 1 and from the ID approximation Equation 6 
for a given 2D rotation profile resulting from a 2D RLS 
inversion of MDI data. The difference is around 12nHz 
for I = 1, 8nHz for I = 2, 6nHz for I = 3. This is 
due to the fact that for low I, Pi {fJ.) extend far from the 
equator and therefore low-degree sectoral modes are more 
sensitive to the latitudinal dependency of the rotation in 
the convection zone than other sectoral modes. 



1. As already mentioned in Section |2.1.| , the rela- 
tion Equation ^ is not valid for low-degree I and 
therefore may not be suited for the determina- 
tion of the core rotation. This may be corrected 
by assuming that the latitudinal dependence of 
the rotation rate (i.e. fij(r) j — 3,5..) is known 
(taken from some previous 2D inversions for ex- 
ample). With this assumption, we can correct 
the observed sectoral splitting prior to inversion 
by adding for each mode the difference betweeri 
the sectoral splittings computed from Equation hi 
and from Equation^. This difference is plottea 
on Figure ^ as a function of the degree / . 

2. MDI data do not provide the sectoral splittings 
for high I but only upon 18 odd indexed a- 
coefficients. This number of coefficients is how- 
ever high enough so that taking their sum as sec- 
toral splittings is a good approximation. Never- 
theless, the error on these sums (called 'trun- 
cated sectoral splittings' in the following) is al- 
ways higher than the error on oi alone. 

Another possibility is to consider the sectoral split- 
tings Z = 2, 3 as oi coefficients at first approxima- 
tion. In this approach we do not need to correct the 
data because Equation ^ is valid even for low de- 
grees. Nevertheless, we can also use other dataset 
(MDI for example) or our knowledge of the latitudi- 
nal dependence of the rotation (taken from a previous 
2D inversion for example) in order to estimate for 
modes I — 2 and 03, 05 for modes 1 = 5. 

With the same approximation as in Equation ^ we 
can write for the a-coefficients; 



(9) 



so that, for both approaches, the solution obtained 
at ro can ever be seen as an average of the rotation 



of the form: 

Sl(ro)= J J[J2 (^o)^ni {r)Wi (/i)) n{r, ii)drdti 

nl 

(10) 

The 2D averaging kernels fdefined by the term in 
parentheses in Equation nO) can therefore be esti- 
mated from the ID averaging kernels Equation |^ ob- 
tained at a target location tq by adding the angular 
part Wi{fj,) of the 2D rotational kernel that corre- 
sponds to the data really inverted. 

Figure shows on the left the 2D averaging ker- 
nels f denned by the term in parentheses in Equa- 
tion [lO) obtained at ro = O.2i?0 by inverting GOLF 



sectoral splittings together with MDI 'truncated sec- 
toral' splittings i.e.: 



for I < 3 



^ Q,,(/i) forZ>3 



(11) 



whereas, on the right, it shows the 2D averaging ker- 
nel obtained at the same target location by invert- 
ing MDI ai coefficients together with GOLF sectoral 
splittings for Z = 1, 2, 3 i.e.: 



for / < 3 



Qzi(/i) = |(l-Ai2) forZ>3 



(12) 



The two corresponding ID averaging kernels (cf. Fig- 
ures. I and |) are the integral over latitude of these 
2D kernels and are very similar: well localized near 
Q.2Rq and without contributions near the surface. 
The angular part of the averaging kernel for sectoral 
splittings strongly depend on the degree I and this 
leads to the oscillatory behaviour of the 2D kernel 
with very high peaks near the surface. From this 
plot it is clear that the interpretation of the result 
obtained by inverting sectoral splittings is possible 
only if we have already a good knowledge of the lati- 
tudinal dependence of the rotation. Furthermore, as 
already pointed out, this knowledge is needed in or- 
der to correct the low-degree sectoral splittings for 
which the ID approximation is not valid. At the op- 
posite, for fli coefficients Q;i(m) is independent of I 
(cf. Equation 12). Therefore the surface oscillatory 
behaviour on tne right panel of Figure ^ comes only 
from the use of sectoral splittings for Z = 2, 3. In this 
case, the knowledge of the rotation profile is needed 
only near the surface in order to interpret the result. 

We can summarize the results of this study in few 
points: 



1. It is clear that it is better to use ai coefficients 
than sectoral splittings when we have them. 

2. In the case of GOLF data, we have access only to 
sectoral splittings for Z = 2, 3. If we want to use 
them, it seems more reasonable to try to correct 
them by doing some assumptions on 03 and 05 



Figure 3. 2D averaging kernels for ID inversions computed at Tq = 0.2Rq (see Section 4)- Left panel: inversion of 
GOLF sectoral splittings for I < 3 with MDI 'truncated sectoral'(see text) splittings for I > 3. Right panel: inversion 
of GOLF sectoral splittings for I < 3 with MDI ai coefficients. The peaks near the surface are truncated on the plot for 
clarity. 



for these modes and to use the exact ID integral 
rather than correcting all the modes in order to 
use the approximated ID integral Equation |4| 

3. In both approaches we can in principle obtain a 
result easy to interpret if the latitudinal variation 
of the rotation is assumed to be known exactly. 
But by inverting ai for I > 3 together with sec- 
toral splittings for ^ < 3, we just have to make 
assumptions on the surface rotation. 

4. Several ways can be followed for correcting the 
results obtained by using sectoral splittings in ai 
inversions. We can either 

• take a guess rotation and integrating its sur- 
face part with the 2D kernel in order to cor- 
rect the solution after the inversion or 

• take 03 and 05 from other datasets or 

• calculate these coefRcients from the guess 
rotation and subtract them from sectoral 
splittings before the inversion. 

The best is probably to compare the effects of 
all these corrections and to compare with the in- 
version of the corrected 'truncated sectoral' split- 
tings. In any case assumptions are needed on the 
latitudinal dependence of the rotation. As this 
dependence can not be known exactly it should 
be interesting to study in future works how these 
assumptions increase the uncertainties on the so- 
lution. 



The next Section shows some preliminary results ob- 
tained with these different approaches for the use of 
the combined MDI and GOLF data. 



5. RESULTS AND DISCUSSIONS 

G0LF(I<3) + MDI144(I> = 3) SOLA (□! coefficients) 
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Figure 4- SOLA inversion of MDI ai coefficients and 
GOLF 1 = 1,2 sectoral splittings. The corresponding ID 
averaging kernels are shown below each point of the solu- 
tion. Vertical and horizontal error bars represent respec- 
tively the error and the resolution defined as the FWHM 
of the averaging kernel. 

In this preliminary work, we have chosen, as a first 
step, to use only / = 1, 2 GOLF sectoral splittings in 
order to reduce the difficulties discussed in the pre- 
vious sections. Figures ^ and |5| show the results ob- 
tained by inverting the ai coemcients together with 
I = 2 sectoral splittings by using SOLA and MOLA 
inverse methods. In SOLA method the trade-off pa- 
rameter is rescaled at each target location to ob- 
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G0NG(I<3) + MDI144(I>=3) UOLA (a1 coefficients) + RLS 




G0NG(i<3) + MDil44(l> = 3) MOW. (al coefficients) + RLS 







];.tVtIII 1 mil 1 jTr^mi n' " ■ - ■ ■ -r ■ ■ 


1 1 1 1 1 1 1 


1 





0.4 0.6 

f/Ro 



Figure 5. MOLA inversion of the same data as in Fig- 
ure 4- Five regularizing parameters have been used at each 
target location. The full line shows an RLS inversion of 
the corrected 'truncated sectoral splittings' (see text) and 
the dotted line gives Isigma errors on this solution. 



G0NG(l<3) + MDI144(i> = 3) MOLA (connected sectoncl splittings) + RLS 




Figure 6. MOLA and RLS inversion of the corrected 
'truncated sectoral splittings '. Low regularization has been 
used in MOLA inversion. 



tain more localized averaging kernels. For MOLA 
method five trade-off parameters fiave been used at 
each, target location. In the case of low regularization 
MOLA averaging kernels have some small oscillatory 
parts. The two methods fall in good agreement down 
to O.2Oi?0 where the rotation rate increases up to 
500 ± 60nHz. Below O.15i?0 both methods fail to 
peak kernels. On Figure we have shown a RLS so- 
lution obtained by inverting sectoral splittings. As 

expected the two solutions (fieq and f2i) differ in the 
convection zone where the rotation rate vary with lat- 
itude. In the core, it is difficult to use RLS method 
because it is a global method and if one try to obtain 
a well localized averaging kernel down to 0.2Rq then 
we have to decrease the regularization and the so- 
lut ion becomes very oscillating evervwhere with big 



Figure 7. The same as Figure 5 but using GONG ai co- 
efficients for 1—1,2,3 instead of GOLF sectoral splittings. 



shown on the plot) computed below this point are 
still localized near O.4i?0 so that, with this method, 
there is no conclusion on the core rotation in terms 
of weighted average of the true rotation. Neverthe- 
less, we can use this solution and look at the residuals 
for each mode. The global normalized of the in- 
version is 1.2. Now, if we look only at low-degree 
GOLF modes, the 'partial normalized x^' is around 
0.5 showing that within error bars GOLF data are 
in good agreement with a constant rotation below 
O.4i?0 and that GOLF errors are probably not un- 
derestimated. Furthermore, looking at the residuals 
for each individual splittings can help in the signal 
analysis by pointing out some modes with high resid- 
uals that may be reanalyzed in order to become more 
confident on the result. 

Following the discussion of the previous sections we 
have also inverted the sectoral (or 'truncated sec- 
toral') splittings corrected by using the latitudinal 
dependence of the rotation found by a 2D RLS inver- 
sion of MDI data. The result is shown on Figure ^ 
in the case of low regularization. As expected the so- 
lution corresponds to the equatorial rotation profile 
as found by the RLS method in the convection zone 
and the error bars increase compared to the use of 
ai coefficients alone. The solution in the core is a 
little bit higher than found by ai inversion but re- 
mains compatible within error bars showing also an 
increasing rotation rate below O.O25i?0. 

We have also tried to include GOLF / — 3 sectoral 
splittings in our inversions. This leads to a more 
important increase of the rotation rate below O.2i?0 
(around 700 ± 100 nHz). But in this case more work 
is needed in order to become more confident in our 
result. In particular, in that case, we have to test 
the effects of the various corrections suggested in Sec- 
tion^ and to look at their influence on the estimation 
of the uncertainties on the core rotation. 

Finally, we have d one the same analvsis with GONG 



err ar bars. With an optimal L-curve choice of the reg- 
ularizmg parameter the solution is constant (445nii z 
see Figure ||) below O.4i?0 but averaging kernels (not 



low-d egree data ( Rabello-Soares fc Appourchaux 
1998|). In the case of GONG data we have ai coetti- 



cients so that, as quoted before, it may be more suited 
to carry an inversion of these coefficients using Equa- 
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tion |6| rather than using the 'truncated sectoral' spht- 
tings. In order to compare the resuh with the pre- 
vious ones, the GONG oi coefficients for low-degrees 
have been used together with the MDI ai coefficients 
of higher degree modes. Figure |7| shows that these 
data tend to produce a slightly decreasing rotation 
rate below Q.3Rq. Therefore there is still a relatively 
important difference between the solutions obtained 
with the different low-degrees data. These differences 
are significant only if the error bars obtained on the 
solutions are not underestimated and may be related 
to the important dispersion of individual splittings 
measurements (cf. Figure |l]). Furthermore, we must 
notice that whereas MDI and GOLF splittings are for 
the same year of observations(5/96-5/97), the GONG 
data are for the year before (5/95-5/96) and therefore 
the results may not be directly compared. There- 
fore this result needs to be confirmed in future works 
and we have also to inverse GONG ai coefficients for 
all the modes which should be a more self-consistent 
dataset. 



6. FUTURE WORKS AND SOME QUESTIONS 



The problem of inferring the core physics remains 
one of the most important still open question that 
can be addressed by helioseismology. It is therefore 
very important, as a first step, to be sure that dif- 
ferent 'inverters' using different inverse methods can 
reach similar conclusions when they use the same 
datasets. This was the goal of this work in collabora- 
tion between 'inverters' within the GOLF team and 
we have shown that the different approaches of the 
inversion are coherent. After this preliminary work 
several questions need to be addressed, for example: 
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• Can we obtain reliable results below O.2i?0 with 
actual datasets? 

• How the results in the core are sensitive to the 
data used for medium I splittings? 

• Can we explain the differences between the re- 
sults obtained with different data by a system- 
atic bias in some splitting measurements or by 
underestimated errors in a few number of indi- 
vidual splittings? 

• Can these results (obtained with inverse meth- 
ods) be confirmed by using forward me 
ten used for the core rotation problem (|Charbon-| 



neau et al. 199^ ; |Gizon 1998]) 



? 
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